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In this paper we present a uniform way to derive families of maps from the corresponding differential equations describing 
systems which experience periodic kicks. The families depend on a single parameter - the order of a differential equation a > 0. 
We investigate general properties of such families and how they vary with the increase in a which represents increase in the space 
dimension and the memory of a system (increase in the weights of the earlier states). To demonstrate general properties of the 
cf-families we use examples from physics (Standard cr-family of maps) and population biology (Logistic a-family of maps). We 
show that with the increase in a systems demonstrate more complex and chaotic behavior. 
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1. Introduction 

Since the first fractional maps (FM) were derived from frac- 
tional differential equations (FDE) |Q]|2]|3][4]], their investiga- 
tion revealed new unusual properties of the fractional dynami- 
cal systems 16] 13 . Probably the most unusual feature of 
the investigated fractional maps is the existence of the new type 
of attractors - cascade of bifurcation type trajectories (CBTT). 

A cascade of bifurcations, when with the change in the 
value of a parameter a system undergoes a sequence of period- 
doubling bifurcations, is a well known pathway of the tran- 
sition form order to chaos (see for example [9|). It is as- 
sociated with the Feigenbaum scaling, Feigenbaum's func- 
tional equation, and period-doubling renormalization operator 
iflOl [TTl [121 [131 [T4l . The scaling property of dynamical sys- 
tems is reflected in two Feigenbaum constants, which are uni- 
versal for large classes of systems and can be computed ex- 
plicitly using functional group renormalization theory. In the 
case of the CBTT the period doubling occurs without a change 
in a system parameter and is an internal property of a system. 
The FMs are much more complicated than the corresponding 
integer maps - they are maps with memory in which current 
values of the map variables depend on all their previous values 

02J[T6l[m[m[lll|20). The first derived from FDEs FMs inves- 
tigated were two-dimensional Fractional Standard Maps (FSM) 
obtained from the fractional Universal Map derived from the 
corresponding FDE. CBTTs were found in all of them but un- 
derstanding of their origin and the necessary and sufficient con- 
ditions of their existence is far from being complet. The fact 
that memory in fractional maps decays slowly (as a power law) 
does not allow to use a "short memory" principle [21 1 and com- 
plicates the investigation of the CBTTs. This is why the devel- 
opment of a simple one-dimensional model with CBTTs is very 
important. 

The simplest integer one-dimensional map which has been 
used as a playground for the investigation of the cascades of bi- 



furcations in regular dynamics is the ubiquitous Logistic Map 
|22|. There were a few attempts to introduce a Logistic Map 
with memory. Bifurcation diagrams in the Logistic Maps with 
exponentially or geometrically decaying memory were consid- 
ered in [18 23 1 and memory related to a numerical integration 
of the FDE was considered in 11201 . The major conclusion is 
that memory increases stability. The maps were not derived 
from FDEs and did not demonstrate CBTTs. It is impossible 
to derive a fractional Logistic Map (FLM) from the equation of 
the fractional Universal Map introduced by Tarasov (for a re- 
view see Chapter 18 from Q) or in a way similar to the way 
in which the Universal Map obtained by considering a system 
with a periodic sequence of (5-function type pulses (kicks). In 
the following we derive the generalized Universal Map (we'll 
call this map the Universal Map omitting for the sake of brevity 
the word generalized and hope it won't cause any confusion) 
by considering kicked systems with time delays. This allows us 
also to derive the FLM in a way similar to the way in which the 
regular Logistic Map can be derived. These maps turns out to 
be the simplest so far one-dimensional maps with the CBTTs. 

2. Regular (Integer) Maps 

2.1. Universal Map (Two-Dimensional) and Standard Map 
The Universal Map 



Pn+l = Pn ~ TG(x n ), X n+ i = x n + Tp„ 



(1) 



can be derived in a way similar to OEJID. Let's consider dif- 
ferential equation 



x + G(x) J] §(L-(n+sj) = 



with the initial conditions: 

x(0) = xo, P (0) = i(0) = po, 



(2) 



(3) 
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were < s < 1. This equation is equivalent to the Volterra 
integral equation of second kind: 

dTG(x(r)) J] <5(--(n + e))(f-T), (4) 

/7=-CO 

which for (n + s)T < t < (n + 1 + s)T has a solution 

n 

x(t) = x + pot-T^ G(x(Tk + Ts))(t -Tk- Ts), (5) 

k=0 

a 

pit) = x(t) = Po -tY j G(x(Tk + Ts)). (6) 

k=0 

With the definitions 

x„ = x(Tn), p n = p(Tn) (7) 
this gives for t=(n+l)T 

n 

x n+l = x Q +po(n+l)T -T 2 YjG(x(Tk+Ts))(n-k+l-s), (8) 



x n+l — x n + Pn+\T, 



(17) 



k=0 



Pn+i = Po ~ T ^ G(x(Tk + Ts)). 



(9) 



fc=0 



Taking the limit e — > and taking into account continuity of 
x(t) (i(f) = p(t) is finite), we obtain 



Xn+l — x 



+ po(n + l)T-T 2 ^G(x k )(n-k+l), (10) 



k=0 
n 

Pn+1 = Po - T 2^ G(x k ), 

k=0 



which can be written in a symmetric form as a map with full 
memory, where all previous states have equal weights in the 
definition of the new state (see for example |20|): 



x n +\ = x + p k , 

k=l 

Pn+1 = Po ~ T^G(x k ). 

k=0 



(12) 



(13) 



Any system with full memory can be presented as a system 
with one step memory (which sometimes is defined as a system 
with no memory). The simplest one-dimensional map with full 
memory 



x n +\ 



(14) 



k=0 



can be written as 

Xn+l = X n + f{x„). (15) 

Similarly, equations ( p"2} and ( [T3] l can be written as a simple 
iterative area preserving (d(p n+ i,x n+ i)/d(p„,x n ) = 1) process 
with one step memory which is called the Universal Map: 



which is, in essence, the relationship between the values of the 
physical variables on the left sides of the consecutive kicks. 

The well-known Standard Map is the particular form of the 
Universal Map with G(x) = K sin x: 



Pn+i = Pn ~ TKsinx, x n+x = x„ + p n +iT. 
One-dimensional Logistic Map 

Xn+l — Kx n (l X n ) 

may formally be written in the 2D form 

Pn+l — ~G[(x n ), X n+ \ — X n + p n +l 

with 

G[(x) — x — Kx(l — x). 



(18) 



(19) 



(20) 



(21) 



It can't be derived from Eqs. ( 16 1 and ( 17 i but may be written as 
the particular form of the Universal Map with two step memory 
(n > 0) 

Pn+l — Pn ~ G(x n ,X n -\), 



X n + 1 — X n + Pn+l, 



where 



G(x,y) = G l (x)-G l (y) 
and the initial conditions 

xo = x-i - Gi(x-i), po = -Gi(x-i) 



(22) 
(23) 

(24) 
(25) 



(11) with an arbitrary x- 



Map ( 22 1 - ( 25 i may formally be derived from the differential 
equation 



x + G(x(t), x(t - 1)) Yj S(f-n-e) = 



(26) 



with the initial conditions (25 1. But this equation at present has 
no physical or biological justification and won't be considered 
in this paper. 

2.2. Logistic Map and Universal One-Dimensional Map 

If we try to use Eq. Q with the first derivative instead of the 
second to derive the ID analog of the Universal Map, we will 
encounter the following difficulty: in the ID analog of the Eq. 
|5]) x(Tk + Ts) will be undefined because Tk + Ts is a point of 
discontinuity for the function x(t). To overcome this difficulty 
let's consider the equation with the time delay 



oo 

x + G(x(t - AT)) J] 5(1 - (n + £)) = 



(27) 



with the initial condition: 



Pn+l = Pn ~ TG(x n ), 



(16) 



x(0) = xq, 



(28) 



where < s < 1 and < A < 1. Then ID analog of Eq. Q 
becomes 



x(t) = Xo -tY j G(x[T(k + s- A)]). 



(29) 



/>=(> 



x[T(k+s-A)] = x(Tk) because x = for t e (T(k+s- 1), 7\jfc+ 
e)). As a result, the corresponding ID map with full memory 
becomes 



X n +\ - X(j 



Tj]G(x k ), 



(30) 



k=0 



which is equivalent to the ID form of the Universal Map with 
one step memory 



x n+\ — x n TG(x n ). 



Logistic Map can be obtained from Eq. ( 3 1 1 by taking 



G(x) 



1 



-\x-Kx(\ -x)]. 



(31) 



(32) 



Differential equation (27i with no time delay, no delta func- 



tions, and G(x) defined by ( 32 1 is one of the most general mod- 
els in population biology and epidemiology E4l l25l . Three 
terms in G(x) represent the growth rate proportional to the cur- 
rent population, the restrictions due to the limited resources, 
and the death rate. Logistic Map appears and plays an impor- 
tant role not only in the life relevant sciences but also in eco- 
nomics, condensed matter physics, and some other areas of sci- 
ence [25 26]. In many cases time delays play an important role 
(see for example Ch. 3 from |25| and Ch. 3 from ll26l ). In the 
simplest case it can be related to the time of the development 
of an infection in a body until a person becomes infectious, or 
to the time of the development of an embryo. The delta func- 
tions model changes that occur as periodically following dis- 
crete events. 

Finally, let us notice that the absence of a time delay is not 
essential in |2]i considered over any finite time interval and both 
the ID and 2D Universal Maps can be derived from the equation 

-£+G(x(t~AT)) <5(--(«+e)) = 0, s > A > 0, i e {1,2} 

n=—oo 

(33) 

in the limit e — > 0: the 2D Universal Map ( [To) , (17i corre- 
sponds to i — 2 and the ID Universal Map d3jj corresponds to 
i=l. 



where e— > 0, < iV - 1 < a < N, and the initial conditions 

( X) D°- k x)(0+) = c k , k= 1, ...,N. (35) 

The left-sided Riemann-Liouville fractional derivative oDfx(t) 
defined for t > ETJ W\ US) as 

D?x(t) = D" I>r a x(t) = 



1 d" f" x(r)dT 
r(n - a) dt" J (t - T) a ~" + 



(n - 1 < a < n), (36) 



D" = d n /dt n , and qI? is a fractional integral. 

This problem for a wide class of functions G(x) can be re- 
duced \4, 28 29, 30 1 to the Volterra integral equation of second 
kind (f > 0) 

x(t) = f C J f-* 

1 f ' J G(x(t - AT)) 



r(a) 



Jo (t — t) ° ') 

Case A = 0, e > 0, a > 1 has been considered by Tarasov in 
El [3] |31 and its consideration consistent with initial conditions 
leads to the equation (f > 0) 



Zc k 
, =i Flo - A - 1 ) 



Ck f- k - 



—L- V ^ + s >\» ®(t-(k + s)T), (38) 
where ®(f) is the Heaviside step function and the map 



x n +i 



2j T(a 



k=\ 



T{a-k+ 1) 



[T(n + l)] a 



-k 



With the introduction 



and 



p (s) (t) = D s t p(t), s = 0,l,-N-2 



(39) 

(40) 
(41) 



it also leads to 



3. Fractional Maps 



Let's consider ( 33 1 in which the derivative of an integer order 
is substituted by a fractional order a derivative. 

3.1. Riemann-Liouville Universal Map 

In the case of the Riemann-Liouville fractional derivative we 
have 

CO 

D?x(t)+G(x(t-AT)) Yj s[--(n+e)) = 0, e > A > 0, (34) 



P {s \t) - X 

k=l 



[t/T-e] 



c k 



(N-s-l-k)\ 



:t N - 



-l-k 



(N-s-2)\ 



Y G(x((k+s)T))(t-kTf- s - 2 , s = 0, 1, .JV-2. 

(42) 



k=0 



and 



N-s-l 
k=l 



Ck 



(N-s-\-k)\ 



[Tin + 1)] 



N-s-l-k 



3 



nJV-J-1 



(N- i-2): 



J] G(**)(n - Jt + l)^" 5 - 



If momentum is defined in a usual way 

= D]x(t), p s (t) = D s t p(t), s = 0, 1...JV - 2, 

then 



(43) 



(44) 



S(f) = Y\ FT" 



J tO'-l-J-fe 



r(or - 5 - jfc) 



[f/r-s] 

E 



G(*((£ + s)T)) 



T(a - 1 - s) (t-(k + s)T) 2+s - a ' 



s = 0,l,..JV-2, t + (k + e)T 



(45) 



and 



Pn+l 



V - Q ,. [r( W + i)f 

^ r(a - s - fe) 



-1-s-i 



r(a-i-l) 



^G^Xn-^+lf^ 2 , s = 0,l,.JV-2. (46) 



/;=() 



The advantages of defining maps using (40 1 and (41 1 are: a) 



Ci in ( |43) are defined as the initial values of the momentum and 
its derivatives; and b) the maps can be significantly simplified 
for computations. For example, for 1 < a < 2 the map equa- 
tions can be written as 



Pn+i = Pn ~ TG(x n ), 
x„+\ = — — V Pk+i V x a (n - k + 1), 

no-) f-^ 



k=0 



where 



V k (m) = m a - k - (m - 1)""% 
and for 2 < a < 3 as 



p\+\ = p\ ~ TG(x„), 



p n+ l = 7>„ + p„-T G(X„), 



(47) 
(48) 

(49) 

(50) 
(51) 



x n +\ — 



Po 



^— [r(n+ l)r 2 + ^ J Vi(n-t+ 1). (52) 



r(or-l) 



k=o 



If A + 0, then x(f) can be continued on f e [-A,0]. Bound- 
ness of x(t) at t — requires c w = and x(0) = 0. The 



map equations (39 1 and (43 1 are still valid in the limit s — > 



(e > A > 0). In the case a — 2 Eqs. (39i and (43 1 converge 



to Eqs. (10 1 and (Hi, while for a — 1 Eq. (39i converges to 
Eq. (fJl). 

We may call Eqs. ( [34} and ( [35} with various map generat- 
ing functions (MGF) G(x) Riemann-Liouville Universal Map 
generating equations (RLMGE). RLMGE generate families of 



maps (c-RL-families) defined by Eqs. (39 1 and (43 1 with vari- 
ous a. 



3.2. Standard and Logistic a-RL-families 

Standard and Logistic a-RL-families of maps are two of the 
most important examples. 

The Standard a-RL-family of maps (FSMRL) is generated 
by the generating function G(x) = K sin x. For a + 1 we as- 
sume K > 0. For T = 1 consistent with the finite initial coor- 
dinate xq — (c N = but cjv_i + and the initial value of x is 
still singular) x and p coordinates of the map for a > 1 can be 
written as 



= X rv ,M +\) a - k --^- Y sm(x k )(n-k+\) a -\ 

j-i Y(a - k + 1) T(a) 

(53) 



N- 1 



it=1 



(AT - k - 1) 



(n + 1) 



N-k-\ 



K 



(N-2): 



Y J sin(Xk)(n-k+ 1) 



N-2 



(54) 



For a < 1 x„ = for all values of n. We consider this map on a 
cylinder (x mod 2n). 

With the same assumptions as for the Standard Map above, 
the x and p coordinates of the Logistic a-RL-family of maps 
with G(x) = x — Kx(l - x) for a > 1 can be written as 



JV— 1 



Xn+l 



1 " 

fy-^ yp k - Kxk(i - xk)](n - k+ i)q: " 1 ' (55) 



N-l 



Pn+l = ^] 

k=\ 



Ck (n + lf- k - 1 



(N-k-1) 



K 



(N - 2)! 



t ~ &it(1 ~ Xk) ^ n ~ k+l y 



N-2 



(56) 



As for the Standard Map, for a < 1 the Logistic RL-family 
produces identical 0. 

3.2.1. Standard a-RL-family for 1 < a < 2 

The regular Standard Map (a = 2) is one of the best investi- 
gated 2D maps OT1I321 . Case 1 < a < 2 has been investigated 
in a series of papers J5] [7] [8) . Fig. [T] reflects the universality 
of the properties of the maps belonging to the same family and 
their dependence on a. 

Low K (K < K c \). Fixed point (0,0), which is a sink for 
1 < a < 2 when 



K < K cl (a) = 



2T{a) 



a I 



where 



V al = Y J (-lf +1 V 1 a (k), 



(57) 



(58) 



for a — 2 converges to the elliptic point which is stable for 
K < K cl {2) = 4. 



4 




Figure 1: Bifurcations in the Standard ff-RL-family of maps: (a). Bifurcations 
in the fractional Standard Map for 1 < a < 2. Below K = K c i curve (0, 0) 
fixed point is stable. It becomes unstable at K = K c \ and gives birth to the 
antisymmetric T = 2 sink which is stable for K c \ < K < K c 2- A pair of T = 2 
sinks with x n+ \ = x„ - n, p a+ \ = —p n is stable in the band above K = K c i 
curve. Cascade of bifurcations type trajectories (CBTTs) appear and exist in 
the narrow band which ends at the cusp at the top right corner of the figure. 
(Xc,Pc) is the point at which the Standard Map's (a = 2) T = 2 elliptic points 
with ,v„ + i = x„ - n, p n+ i = —p n become unstable and bifurcate. In the area 
below K c 3 (above the CBTT band) the chaotic attractor is restricted to a band 
which width is less than 2n. On the upper curves and above them the full phase 
space is chaotic; (b). Bifurcation diagram for the one-dimensional Standard 
Map {59). 

The one-dimensional Standard Map (a — 1) can be written 

as 

Xn+i = x n - K$m(x„), (mod 2jt), (59) 

which is a particular form of the Circle Map with zero driving 
phase. It has attracting fixed points 2nn for < K < K c \{\) = 2 
and ji + 2nn when -2 < K < (see Fig.fT]?). 

Antisymmetric period 2 point (K c \ < K < K c z). 1 < a < 2 
antisymmetric period 2 sink 

Pn+i = ~Pn, x n+ \ = -x n (60) 

exists for K > K c \(a) and is stable for 

K cl (a) <K< K c2 (a), where K c2 (a) = (61) 

Val 

For the regular Standard Map it gives a well known result 
K c i(2) = 4 < K < K c2 (2) = 2n. For a = 1 stability analy- 
sis produces K c \(\) — 2 < \K\ < K C 2(2) = n for the stability of 
the period two sink x n+1 = 2nm - x„. 

x n +\ = x n + n period 2 point (K > K c o below the CBTT 
band). For 1 < a < 2 two periodic sinks 

Pn+l = ~Pn, X n+ ] = x„ + n (62) 

appear at K — K c i and exist for K > K C 2 below the CBTT 
band. For the Standard Map the corresponding elliptic points 
are stable when 2ti < K < 6 .59. Fo r a — 1 x n+ \ — x„ + n sinks 
are stable when n < \K\ < y/n 2 +2 a 3.445. 
Cascade of bifurcations band. 

At K =s 6.59 in the Standard Map an elliptic-hyperbolic point 
transition, when T = 2 points become unstable and stable T — 4 
elliptic points, appear occurs. Further increase in K results in 
the period doubling cascade of bifurcations which leads to the 
disappearance of the corresponding islands of stability in the 




Figure 2: Cascade of bifurcations type trajectories: (a), a = 1.65, K = 4.5; one 
intermittent trajectory in phase space; (b). Time dependence of the coordinate 
x (x of n) for the case (a); (c). a = 1.98, K = 6.46; zoom of a small feature 
for a single intermittent trajectory in phase space; (d). a = 1.1, K = 3.5; a 
single trajectory enters the cascade after a few iterations and stays there during 
500000 iterations. 



chaotic sea at K w 6.6344. The cusp in Fig. [ij points to ap- 
proximately this spot (a — 2, K w 6.63). Inside the band, lead- 
ing to the cusp, a new type of attractors appears: cascade of 
bifurcations type trajectories (see Fig.|2|. In the CBTTs period 
doubling cascade of bifurcations occurs on a single trajectory 
without any change in the map parameter. A typical CBTT's 
behavior is similar to the Hamiltonian dynamics in the pres- 
ence of sticky islands: occasionally a trajectory enters a CBTT 
and then leaves it entering the chaotic sea (Fig.|2^,b). As a de- 
creases the relative time trajectories spend in CBTTs increases. 
Near the cusp (Fig. |2j;) CBTTs are barely distinguishable and 
the relative time trajectories spend in CBTTs is small. When a 
is close to one, a trajectory enters a CBTT after a few iterations 
and stays there over the longest computational time we were 
running our codes - 500000 iterations. For a — 1 sequence 
of bifurcations: T = 4 at K * 3.445, T = 8 at K * 3.513, 
T — 16 at K x 3.526, and so on leads to the transition to chaos 
at K « 3.532. Single antisymmetric T = 2 (K = 2.5), T = 4 
(K = 3.5), and two chaotic trajectories (K = 4 and K = 5) are 
presented in Fig. [3] 

K < K c3 above the CBTT band - x-proper attractor 
In the one-dimensional map with K > the full phase space 
becomes involved in chaotic motion (we'll call this case "im- 
proper attractor") when the maximum of the function /k(x) 
x K sin x is equal to n which occurs at K max = 4.603339 when 
Xmax = -1.351817. This is the point at which K = K c j,(a) curve 
on Fig.[lJ hits the a — 1 line. In the area between K = K c j,(a) 
curve and the upper border of the CBTT band the fractional 
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Figure 3: Attractors in one-dimensional Standard Map - x„ vs. x n +\ plots: (a). 
A' = 2.5; antisymmetric T = 2 sink; (b). = 3.5; T = 4 trajectory; (c). A" = 4; 
proper attractor (width of the chaotic area is less than 2n); (d). K = 5; improper 
attractor (width of the chaotic area is 2n). 




3.14 



Figure 4: Proper and improper attractors. 3000 iterations on ten trajectories 
with the initial conditions xq = 0, po = 0.001 + 1.65i, i = 0,1,. ..9: (a). A 
proper attractor for K = 4.2, a = 1.1; (b). An improper attractor for K = 4.4, 
a = 1.1. 



attractors are proper (see Fig. |4]i 

K = 2m accelerator mode attractor 

The regular Standard Map has a set of bands for K above 
2m of the accelerator mode sticky islands in which momen- 
tum increases proportionally to the number of iterations n and 
coordinate increases as n 2 . The role of the accelerator mode 
islands (for K above 2n) in the anomalous diffusion and the 
corresponding fractional kinetics is well investigated (see, for 
example, I331I34T ). 

In the one-dimensional map the corresponding bands demon- 
strate cascades of bifurcations (see Fig.fTJi) for above 2n\n\. 
The acceleration in those bands is zero and x increases propor- 
tionally to n. 

The case 1 < a < 2 is not fully investigated. The acceler- 
ator mode islands evolve into the accelerator mode (ballistic) 
attracting sticky trajectories when a is reduced from 2 for the 
values of K which increase with the decrease in a (Fig. [S]}). 



Figure 5: Accelerator mode attractors. 25000 iterations on a single trajectory 
with the initial conditions xq = 0, po = 0.1: (a). CBTT-type accelerator mode 
attractor for K = 5.7, a = 1.03; (b). Accelerator mode attractor for K = 7.6, 
a = 1.97. 



When the value of a increases from 1, the corresponding ballis- 
tic attractors evolve into the cascade of bifurcation type ballistic 
trajectories (see Fig. [5^) for the values of K which decrease with 
the increase in a. This could mean that corresponding features 
in the one- and two-dimensional maps (at least for K = 2n) are 
not connected by the continued change in a. 

3.2.2. Standard a-RL-family for 2 < a < 3 
a = 3 Universal and Standard Maps 



Following Eqs. (|39]l,(42 1, and (43 1, and assuming T = 1 the 
three-dimensional universal map can be written as (y — p, z — 
P) 

1 " 

x n+ x = Z ^(n+l) 2 +y Q (n+l) + x --Y J G{x k )(n-k+\)\ (63) 

k=0 

y n+ \ = zo(n + 1) + y - ^ G(x k )(n - k + 1), (64) 



k=0 



Zn+l = ZQ ~ ^ G(x k ), 

which can be reduced to 



(65) 



k=0 



or 



y n+ x = -G(x n )+y n +z„, 

Zn+\ = -G(x n ) +Z„, 
1 

y n +\ = y n + Z„-H, 
z„+i = -G(x n ) +z„, 



(66) 



(67) 



which is a volume preserving map. 

Let's assume that this map has a fixed point (xo,yo,Zo)- Then 
Zq — yo — G(xq) = and stability of this point can be analyzed 



-2 





x 



Figure 6: Phase space of the 3D Standard Map {68} for K = 3: (a). Phase 
space in three dimensions; (b). Projection of the phase space on the x-y plane. 



considering the eigenvalues of the matrix (corresponding to the 
tangent map) 



( 1 - 0.5G(* ) 1 0.5 ) 
-G(x ) 1 1 
-G(x ) 1 



The only case in which the fixed point is stable is G(xq) = 0, 
when Xi = A2 = A3 = 1. 

For the three-dimensional Standard Map with G(x) = 
K sin(x) 

x„+\ = x n +y n+ \ - -Zn+i, (mod 2n), 



y n +i =y n +Z«+i. (mod2;r), 
z„+\ = -Ksm(x n ) + z n , (mod An) 



(68) 



it means that fixed points (2nn,2nm,4nk) and (2nn + 
n, 2nm, 4nk) , n e Z, m e Z, k e Z, are unstable for all ^ 0. 
Ballistic points sin(jt) = -Ann, y - 2nm, z = 47r£, which 
appear for > An, are also unstable. 

Stability of T = 2 points is defined by the eigenvalues of the 
matrix 



C 1 -0.5/Tcosjci 1 0.5 



-/Tcos xi 
-A' cos 



1 





( 1 -0.5Kcosx 2 1 0.5 



-K COS ^2 
-/T COS X2 



1 





For the stable (on the torus) period two point 

zu yi = y - n(2n + 1), 7irsinxi=2zi, 

Z2 = -Zi. = -y - ?r(2n + 1), X2 = Xi - n(2n - 1), (69) 
where neZ, the eigenvalues are 

{l, -(8 -/T 2 cos 2 x, i/Tcosxi ^K 2 cos 2 X] - 16)} (70) 



Stable T - 2 points exist along a line defined by Eqs. ( |69| > for 
all values of z satisfying the condition 



K 2 - 16 < 4z 2 < 



(71) 



An example of phase space for K = 3 in three dimensions and 
its projection on the x-y plane is given in Fig. [6] In this case T = 
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Figure 7: FSM for 2 < a < 3: (a). 3D phase space for AT = 1, a = 2.01 
obtained on a single trajectory with xq = po = and = 0.01; (b). Projection 
of the phase space in (a) on the x-y plane; (c). Projection of the phase space for 
K = 0.2, a = 2.01, xq = pq = on the x-y plane obtained using 20 trajectories 
with different initial values of p l Q \ (d). The same as in (c), but for K = 4 and 
a = 2.9. 



2 points are stable for -1.5 < z < 1.5 and the space around the 
line of stability presents a series of islands (invariant curves), 
islands around islands, and separatrix layers. As K goes to zero, 
the volume of the regular motion shrinks. For small values of K 
the line of the stable T — 2 points exists for —K/2 < z < K/2. A 
different form of the 3D volume preserving Standard Map was 
introduced and investigated in details in [35 1. 
2 < a < 3 

Assuming G(x) = Ksin(x) and T = 1 in Eqs. (|50l>-(|52]>, the 
FSMRL for 2 < a < 3 can be written as 

P l n+ i =p 1 n -Ksin(Xn), (72) 
Pn+i = Pl+ Pn-K sin(x„), (mod 2n), (73) 



x n +\ — 



PO 



F(a - 1) 



n + ir 2 + Z ^ y « ( " _ * + 1} ' (74) 



(mod 2n). 

In our simulations we did not find a stable fixed point even for 
small values of K (see Figs. |T| (c). Simulations show that for 
this map there are attractors in the form of the attracting multi- 
period lines with constant x (see Fig. [7] (a), (b), and (d)). For 
most of the values of the map parameters the phase space is 
highly chaotic. 

This case and the transition from the 2D Standard Map to the 
3D Standard Map, as well as 3D Standard Map, is not yet fully 
investigated. 

3.2.3. Logistic a-RL-family for 1 < a < 3 

Case a — 1 is probably the best investigated map - the regular 
Logistic Map. For 1 < a < 2 it can be written as 



which for a 



Pn+l - Pn + Kx n (l - X n ) - X„, 
1 " 

x„+i = 7T— - V Pk+\Vl{n - k + 1), 
- 2 can be reduced to 

Pn+l = Pn + Kx n (l - X n ) - X„, 
x n+l — x n + Pn+l ■ 



(75) 
(76) 

(77) 
(78) 



This map is similar to the original area preserving quadratic 
map considered by Henon ll36l (for a recent review see book 
1 37]). As for the Standard or-RL-family, for the Logistic a-RL- 
family of maps we are interested in the stable fixed and periodic 
points in order to follow their evolution as a function of a. For 
-3 < K < 1 (0, 0) fixed point of the 2D Logistic Map is stable. 
At K = 1 it becomes unstable but fixed point ((K - l)/K,0) 
becomes stable. This point is stable for 1 < K < 5. At K = 
5 a bifurcation occurs and this stable fixed point changes to a 
couple of T — 2 stable points 



K + 3+ y/(K + 3)(K - 5) 
x — , v 



For 2 < a < 3 the map can written as 



Pn+l = Pn + Kx„(l - x n ) - x„, 



Pn+\ = P\ + Pn + Kx„(l - X n ) ■ 



V(T + 3){K - 5) 



2K 



(79) 

(80) 
(81) 



X n +l — 



P0 



T(a-1) 



n + ir- 2 + -^- fjP l k+l V ] a (n-k + l), (82) 



which for a — 3 can be reduced to 



X n +\ — x n +y«+l 2^"+l' 



y n +i = yn +Z„+U 

Zn+l = Kx„(l - X n ) ~X„+ Z n - (83) 

3D quadratic volume preserving maps were investigated in 



|. Everything stated in Subsection 3.2.2 for the 3D Uni- 
versal Map is still valid for the 3D Logistic Map. More on the 
Logistic ff-RL-family of maps will be presented in the part II of 
the article. 

3.3. Caputo Universal Map 

The Logistic and Standard a-RL-families give identical ze- 
ros in the case a < 1 . This is not the case for the Logistic and 
Standard a-Caputo-families (FLMC and FSMC), for which the 
condition xq — can be removed. In order to get some insight 
into the origin of one of the most interesting features of frac- 
tional maps - CBTTs, we complete this part of the paper with 
the analysis of the simplest case of CBTTs in the Logistic and 
Standard a-Caputo-families of the maps when a < 1 . But let's 
first introduce the Caputo Universal Map. 



Similar to ( 34 1 in the case of the Caputo fractional derivative 
we have 

oo 

$Dfx(t) + G(x(t - AT)) ^--(» + «)) = 0, b > A > 0, 

(84) 

where s — > 0, < TV - 1 < a < N, and the initial conditions 

(D?Jt)(0+) = b k , k = 0,...,N- 1. (85) 

The left-sided Caputo fractional derivative {j Dfx(t) defined for 
t> 0112111271 1281 as 



%D?x(t)= I"'" D n t x{t) = 



1 C' D'^ X (T)dT 

r(n - a) Jo (t - T) a -" +1 (n 



1 < a < n). (86) 



This problem is equivalent to the following Volterra integral 
equation (t > 0) 11411281 

k=o K - 

1 C' , G(x(t-AT)) 



a) Jo (f-r) 1 - ^ v r '> 



r( 



Then we also have for t > 



k=a 



T [,/ ^ B] G(x((k + e)T)) 



- Y 

(a) 4-i 



T(a) (t-(k + s)T) 



^®{t-(k + s)T) (88) 



and the map 



|V-I 



k=0 



Xn+i = Yj% [T{n+l)]k 



yd' 



— YG(x k )(n-k+lT-\ 



(89) 



which for < a < 1, assuming T = 1, gives 



1 " 



3.3.7. Standard and Logistic a-Caputo-families for < a < 1 
For < a < 1 the Standard a-Caputo-family can be written 

as 

x„ = xq - — V S1 ° ,^ , (mod 2;r) (91) 
r(a) ^ (n-^) 1 ff 

and the Logistic a-Caputo-family as 



n— I 



1 Y - x k ) - x k 

X n ~ Xq + —— > — — . (VZ) 

r (a) ^ (« - ^) J ff 
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Figure 8: Cascades of bifurcations and CBTTs in Standard and Logistic a- 
Caputo-families (FSMC and FLMC) for < a < 1 : (a). Cascade of bifurca- 
tions for the FSMC J5TJ with a = 0.8 obtained after performing 10000 iter- 
ations on a single trajectory with xq = 0.1 for various values of K; (b). The 
same as in (a) but for a = 0.3; (c). The same as in (a) but for a = 0.01; (d). A 
single CBTT for the FLMC with a = 0.4, K = 6.075, and the initial condition 
xo = 0.1; (e). The same as in (a) but for the FLMC (92) ; (f). The same as in (b) 
but for the FLMC; (g). The same as in (e) but for a = 0.1; (h). The same as in 
(g) but after 100 iterations. 



As one may see, from (91 1, (92 1, and Figs.[8](a)-(c) and (e)-(g) 
a decrease in a and corresponding decrease in the weights of 
the earlier states (decrease in the memory effects) leads to the 
shift of the bifurcation curves to the higher values of the maps' 
parameter K, which means increased stability. When a — > 1 the 
bifurcation curves converge to the curve in Fig.JTJ) for the FSM 
and to the well known bifurcation curve for the Logistic Map in 
the case of the FLM. 

For the values of K at which periodic T > 2 stable points 
exist the individual trajectories are the CBTTs (see Fig. [8}l), 
where an increase in the number of the map iterations leads 
to the change in the map's stability properties. These are the 
simplest so far examples of the maps with the CBTTs. As a 
result, the FSMC and the FLMC can be used as the simplest 
models for studying the CBTT phenomenon. 

The presence of CBTTs leads to the dependence of the the 
bifurcation curve on the number of iterations. From Figs.[8](g) 
and (h) one may see that with the increase in the number of 
iterations in the the FLMC the bifurcation curve becomes more 
pronounced and shifts to the left (to the lower values of the map 
parameter K). 

4. Conclusion 

Differential equations with the integer order derivatives and 
the terms containing periodic delta-function shaped kicks are 
frequently used in science to model the behavior of real sys- 
tems. They can be easily integrated over an integer number of 
periods in order to obtain maps, which then can be numerically 
investigated in order to derive the general properties of the orig- 
inal systems. Well known examples of such maps are the Stan- 
dard Map in physics and the Logistic Map in population biol- 
ogy. In this paper we proposed continuation of the integer maps 
into the area of the fractional maps corresponding to the arbi- 
trary positive values of the orders of the fractional derivatives 
in the original differential equations. We introduced the notions 
of the Riemann-Lioville and Caputo a-families of maps, which 
depend on a single parameter - order of the fractional derivative 
a > 0. Increase in a corresponds to the increase in the dimen- 
sion of the space (which can be important in physical systems) 
and to the increase in the memory effects (the weights of the 
earlier states increase with the increase in a), which is impor- 
tant in population biology (information about past states of the 
system can be recorded in the DNA of individuals, the society 
as a whole can regulate evolution of the population by introduc- 
ing laws based on the past experience). 

Our first results show that increase in a leads to more a more 
complex behavior of the system. As a — > any memory 
of the past disappear and the system stabilizes at zero for the 
Riemann-Lioville families or at constant value for the Caputo 
families. For < a < 1 the behavior of the systems is qualita- 
tively similar to the behavior of the system with a = 1 with sta- 
bility reducing as a —> 1. Within the bifurcation (T > 2) band 
of the parameter values fractional systems evolve as CBTTs. 

For 1 < a < 2 the systems become two-dimensional and their 
phase space contains different kinds of arttractors like sinks, 
chaotic attractors, CBTTs, and intermittent CBTTs. 
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Case 2 < a < 3 requires more thorough investigation. 
Our preliminary results show the existence of very complicated 
three-dimensional attractors and instability of the fixed points. 
In general, as a increases, the systems become more and more 
chaotic. 

In the second part of this paper we'll present more results on 
the fractional Logistic Map with a > 1, the fractional Standard 
Map with a > 2, and the CBTTs. 
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